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Interacting dark matter and modified holographic Ricci dark energy 

induce a relaxed Chaplygin gas 
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We investigate a model of interacting dark matter and dark energy given by a modified holographic 
Ricci cutoff in a spatially flat Friedmann-Robertson- Walker (FRW) space-time. We consider a 
nonlinear interaction consisting of a significant rational function of the total energy density and its 
first derivative homogeneous of degree one and show that the effective one-fluid obeys the equation 
of state of a relaxed Chaplygin gas. So that, the universe is dominated by pressureless dark matter 
at early times and undergoes an accelerated expansion in the far future driven by a strong negative 
pressure. We apply the ^-statistical method to the observational Hubble data and the Union2 
compilation of SNe la for constraining the cosmological parameters and analyze the feasibility of 
the modified holographic Ricci cutoff. By using the new Om diagnostic method, we find that the 
effective model differs substantially from the A-CDM one, because it gets the accelerated expansion 
faster than the A-CDM model. Finally, a new model with a third component decoupled from the 
interacting dark sector is presented for studying bounds on the dark energy at early times. 

PACS numbers: 
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I. INTRODUCTION 

The holographic principle states that the maximum 
number of degrees of freedom in a volume should be 
proportional to the surface area [l[, @], [3], [J]. Using 
the effective quantum field theory it was shown that the 
zero-point energy of a system with size L should no ex- 
ceed the mass of a black hole with the same size, thus 
L 3 p\ < LMp, where p\ corresponds to the quantum 
zero-point energy density [5[ and Mp 2 = SttG. The lat- 
ter relation establishes a link between the ultraviolet cut- 
off, define through p\, and the infrared cutoff which is 
encoded by the scale L. By applying this novel prin- 
ciple within the cosmological context, one can assume 
that the dark energy density of the universe p x takes the 
same form of the vacuum energy, pa — Px- Taking the 
largest L as the one saturating the above inequality, it 
turns out to be the holographic dark energy is given by 
p x = 'ic 2 M 2 p L~ 2 , where c is a numerical factor. Here, 
the infrared cutoff will be considered as a function of the 
cosmic time so the holographic dark energy will evolve 
dynamically and in order to be taken as a good candi- 
date for occupying the place of the cosmic fuel called dark 
energy it should be able to mimic the current accelerated 
expansion phase of the Universe. So far in the literature, 
the IR cutoff has been taken as the large scale of the 
universe, Hubble horizon [fj|, [Z| , particle horizon, event 
horizon or generalized IR cutoff ||, d, OS, & 0- 

An interesting property of the modified Ricci cutoff (lOj 
refers to the fine-tuning problem, that is, due to the holo- 
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graphic dark energy density is based on spacetime scalar 
curvature without involving a Planck or high physical 
energy scale the fine-tuning problem might be avoided, 
further the coincidence problem seems to disapear within 
this context Q. For all these reasons, we will focus on 
a holographic model where the dark energy den sity cor- 
responds to the modified Ricci cutoff [H|, [13], [14[ fl5J |. 
[la ] extending the recently investigated linear interaction 
model [16| to the case of nonlinear ones. 

Another central aspect of the modern cosmology is as- 
sociated with the exchange of energy between the dark 
components, i.e., the dark matter not only can feel the 
presence of the dark energy through a gravitational ex- 
pansion of the Universe but also can interact between 
them. Because the nature of the dark sector still remains 
unknown there is no a fundamental reason to specify the 
coupling between dark energy and dark matter. In fact, 
one follows a phenomenological approach by studying the 
properties hidden in each type of interaction and then one 
confronts the theoretical model with the available obser- 
vational data. Reference [l7| investigated several linear 
and nonlinear interactions in the dark sector, generaliz- 
ing previous works on the literature. It was introduced 
an effective one-fluid description of the dark components 
and shown that interacting and unified models are related 
to each other. Also, a generic nonlinear interaction in- 
duces an effective equation of state resembling a variable 
Chaplygin gas model, giving rise to the "relaxed Chap- 
lygin gas model" [17| ■ Within this framework, we are 
going to explore an holographic model and its relation to 
relaxed Chaplygin gas. 

The paper is outlined as follows. In Sec. II we intro- 
duce the model and discuss about the exchange of energy 
in the dark sector where the modified Ricci cutoff is taken 
into account. For the selected nonlinear interaction term, 



we show that the effective one-fluid obeys the equation 
of state of a relaxed Chaplygin gas. In Sec. Ill we find 
the most probable value of the cosmological parameters 
by applying the x 2 minimization method to the Hubble 
data and the Union2 compilation of SNe la. We use the 
best-fit value of the parameters for confronting our model 
with the standard A-CDM and apply the Om diagnostic 
method. Finally, in Sec. IV we summarize our findings. 



II. THE MODEL 



with A7 = a 
given by 



/3, while the total pressure p = p c + p x is 



P(P,P') 



-P~P 



(7) 



From now on we will use the MCE (UJ instead of the 
WCE with variable u x because it is simpler, and in- 
troduce an interaction between both dark components 
through the term 3HQ M into the MCE (0} with con- 
stant coefficients, so 



p c + ap c — -Qm, 



(8) 



We present an holographic cosmological model, where 
we associate the IR cutoff with the modified Ricci's radius 
and take L~ 2 in the form of a linear combination of H 
and H 2 . After that, the modified holographic Ricci dark 
energy (MHRDE) 0, p x = 3c 2 M 2 P L~ 2 , becomes 



a — (3 



H 



3a 



-H' 



(D 



Here H = a/ a is the Hubble expansion rate, a is the scale 
factor and a, (3 are free constants. Introducing the vari- 
able 77 = ln(a/ao) 3 , with ao the present value of the scale 
factor and ' = d/drj, the above MHRDE (fTJ) becomes a 
modified conservation equation (MCE) for the cold dark 
matter p c and the MHRDE 



p' = -ap c - /3p x , 
after using the Friedmann equation, 
3H 2 = p c + p x , 



(2) 



(3) 



for a spatially flat FRW cosmology and p — p c + p x . The 
MCE © looks as it were a conservation equation for 
both dark components with constant equations of state. 
In connection with observations on the large scale struc- 
tures, which seems to indicate that the Universe must 
have been dominated by nearly pressurclcss components, 
we assume that p c includes all these components and has 
an equation of state p c = while the MHRDE has a 
barotropic index uj x = p x / ' p x , so that the whole conser- 
vation equation (WCE) becomes 



Pc 



l)Px- 



(4) 



The compatibility between the MCE © and the WCE 
(d|) yields a linear dependence of the equation of state of 
the MHRDE 



0J X 



(o-l)r + G8-l), 



(5) 



with the ratio of both dark components r = p c /p x - 
Solving the linear algebraic system of equations @ and 
p = p c + p x we obtain both dark energy densities as func- 
tions of p and p' 



Pp + P' 

A7 



Px 



ap 



A7 



(6) 



Px + PPx = Qm- 



(9) 



Finally, from Eqs. (JS]) and ©, we obtain the source 
equation [17[ for the energy density p 



p" + (a + P)p' + a/3p = QAf A7. 



(10) 



Now we consider cosmological models where the inter- 
action Qm between both dark components is nonlinear 
and includes a set of terms which are homogeneous of de- 
gree 1 in the total energy density and its first derivative 



Qm = 



(a/3-1) (a + J3 



2) 



A7 



A7 



vp 

pA 7 ' 



(11) 



where v is a positive constant that parameterizes the in- 
teraction term Qm- Replacing (fTTj) into ([TU| it turns into 
a nonlinear second order differential equation for the en- 
ergy density 



pp" + (2 + v)pp' + vp 12 + p 2 = 0. 



(12) 



Introducing the new variable y = p^ l+v > into the latter 
equation one gets a second order linear differential equa- 
tion y" + (2 + v)y' + (1 + v)y = 0, whose solutions allow 
us to write the energy density 



P 



Pwa 



p-2oa 



-3(n 



1/(H 



(13) 



where pio and P20 are positive constants. From Eqs. (|6]) 
and ([7]) we have both dark energy densities and the total 
pressure 



a-P 



13-1 + 



{l + u){l + p2oa- Zv /p 10 ) 



a-P 



{l + v){l + P2 a- 3l/ /pw) 



V 



1 



vp w _ a_ 
v p v 



(14) 



(15) 
(16) 



From these equations we see that an initial model of in- 
teracting dark matter and dark energy can be associ- 
ated with an effective one-fluid description of an uni- 
fied cosmological scenario where the effective one-fluid, 



with energy density p = p c + p x and pressure (|16|) . 
obeys the equation of state of a relaxed Chaplygin gas 
p = bp + f(a)/p u : where b is a constant [13] ■ The effec- 
tive barotropic index uj = p/p = uj x p x / p reads, 



vp w 



{l + v){p w +p 20 a-* v )' 



(17) 



At early times and for v > 0, the effective energy den- 



sity ([13]) behaves as p 



,-3 



, the effective barotropic 



index (fT7| 7 » 1 and the effective fluid describes an 
Universe dominated by nearly pressureless dark matter. 
However, a late time accelerated Universe i.e., uj < —1/3 
with positive dark energy densities require that v > 1/2, 
/? < 1 and a > 1. From now on we adopt the latter 
restrictions. 

In the case v = 1, one gets an equation of state (fl6|) 
similar to a variable Chaplygin gas as the one studied 
in [l3|, [19|. In these works the authors have used the 
gold sample of type la supernova data as well as the x- 
ray gas mass fractions in galaxy clusters to obtain some 
observational constraints of the model. 



III. OBSERVATIONAL HUBBLE DATA 
CONSTRAINTS 



Now, we are going to find the observational constraints 
on the parameter space of the MHRDE using the observa- 
tional H(z) data. The function H(z) plays a crucial role 
in understanding the properties of the dark energy, since 
its value is directly obtained from the cosmological ob- 
servations. Using the absolute ages of passively evolving 
galaxies observed at different redshifts, one obtains the 
differential ages dz/dt and the function H(z) can be mea- 
sured through the relation H(z) = —(1 + z)~ 1 dz/dt. We 
take the data coming from the Gemini Deep Deep Survey 
[Ilarchival data (^M, pg] , US, (13|, HE|, H3] , Simon et 
al [23] and compare it with the H(z) to attempt to obtain 
the constraints for the current cosmological models. The 
12 observational H(z) data from [28[ are listed in Table 
I. There, H b s (zi) and H b s (zk) are uncorrelated because 
they are obtained from the observations of galaxies at 
different redshifts. 

The probability distribution for the ^-parameters is 
given by a Gaussian density (see e.g. [23]) 



P(0) = Ze 



-X 2 {6)/2 



fc=l 



[H(e,z k )-H obs (z k )}' 
a(z k f 



where Z is a normalization constant, H b s (zk) is the ob- 
servational H(z) data at the redshift Zk, cr(zk) is the cor- 
responding la uncertainty, and the summation is over the 
12 observational H(z) data of Table I. Since we are inter- 
ested in obtaining the bounds for the model parameters, 



we have adopted as a prior H$ — 72. 2 ±3. 6 km s l Mpc l 
[301 ] . The Hubble expansion of the model becomes: 



H(0; z) = H {B(1 + zf + (1 - B)(l + Z )*»+V}*™ 

(19) 



B[9) = 



1 



[ a (n x0 - 1) + (1 - /3n x0 )} (20) 



where 9 = {a, /?, fl x o, v) and we have used that P02/ poi ~ 
(B — l)/B. The two independent parameters a and /? will 
be fixed along the statistic analysis. Then, for a given 
pair of (af,[3f), we are going to perform the statistic 
analysis by minimizing the x 2 function to obtain the best 
fit values of the random variables 8 C = {v, Q x q} that 
correspond to a maximum of Eq. (fl~8| . More precisely, the 
best-fit parameters 9 C are those values where Xmin(®c) 
leads to the local minimum of the x 2 (#) distribution. If 
Xdof = xl lin ( e c)/(N - n) < 1 the fit is good and the 
data are consistent with the considered model H(z;9). 
Here, N is the number of data and n is the number of 
parameters [29|]. The variable \ 2 ' IS a random variable 
that depends on TV and its probability distribution is a 
X 2 distribution for N — n degrees of freedom. 

Besides, 68.3% confidence contours in the (1/, f2 x o) 
plane are made of the random data sets that satisfy the 



inequality Ax = X (#) 



A ; 



: ) < 2.30. The lat- 



ter equation defines a bounded region by a closed area 
around 9 C in the two-dimensional parameter plane, thus 
the la error bar can be identified with the distance from 
the 9 C point to the boundary of the two-dimensional 
parameter plane. It can be shown that 95.4% confi- 
dence contours with a 2a error bar in the samples sat- 
isfy Ax 2 < 6.17 while the data within 99.73% confidence 
contours with a 3er error bar are accommodated in the 
domain defined by Ax 2 < 11.8. After performing this 
analysis we arc in position to get confidence contours 
in the (1/, £l x o) plane, thus using the x 2 (a/,/3/, v, Q, x q) 
distribution one can find the 68.3%, 95.4%, and 99.73% 
confidence contours respectively [see Fig.(fT|)]. We have 
taken the point of reference (a/,/3/) = (1.01,0.15) but 
it is possible to show a wide set of admissible values 
for a and f3 which leads to a good fit [see Table fllf ]. 
Thus, from Fig.([T]) we get the best fit at 9 C = {y, O. x o) = 
(1.19 ± 0.12; 0.61 ± 0.02). It corresponds to a local mini- 
mum Xmin = 7-86 leading to a good fit with x\ f = 0.786 
per degree of freedom. The Ricci's cutoff (a, (3) = (4/3,1) 
does not guarantee the convergence of the minimization 
process. However, the values of {y, £l x o) obtained from 
an holographic dark energy p x oc R. namely (4/3, /3), 
(18) fulfills the goodness condition Xdof < 1- Th° vamcs 
of fl X Q, which varies from, 0.58 to 0.69, do not devi- 
ate significantly from the observational limits provided 
by the WMAP-7 project [H[ with n x0 = 0.73 [see Ta- 
ble (UJ)]. Comparing the Ricci model with the one aris- 
ing from MHRDE for (a = 1.01,/? = 0.15), the former 
gives (v, Oko) = (1.19,0.69), whereas the latter yields 
(is,n x o) = (1.19,0.61), so the Ricci model seems to be 



statistically favored by H(z) data showing a £I x q closer 
to the observational bound reported by the WMAP-7 
project [311 ]. 

We estimate the best value of Ho and fl X Q for the 
A-CDM model using the Hubble data as well as the 
Union2 data for SNe la 32[ . The former dataset leads to 
H = 73. 60 ±3. 18 km s" 1 Mpc" 1 and Q x0 = 0.730 ±0.04 
with Xdof = 0.770 whereas the latter one gives H 



0.73 along with Xdof 



of 

70 kms-'Mpc" 1 and fi x0 - u .,u a M . 6 wiun Xdof 
0.978. Now, in odcr to make possible a comparison with 
the A-CDM model we need to estimate the same types 
of parameters so we take as priors v = 1.19, a = 1.01, 
and /3 = 0.15 but we allow Hq and fl x o as free param- 
eters to be found under the minimization process. It 
turns out that H = 73.011 ± 2.97 km s^ 1 Mpc^ 1 and 
fl x0 = 0.617 ± 0.001 with x\ of = 0.779 < 1 [sec Fig.©], 
then both models give cosmological bounds of the pai r 
(Ho, Qox) very consistent with the those reported in [31| . 
In order to compare the Hubble data (12 points) with 
the Union2 compilation of 557 SNe-Ia [32[ we proceed as 
follows; thus, we took as priors Ho = 72.2 km s~ J Mpc -1 , 
a = 1.01 and (3 = 0.15 in both cosmological data. We 
found the best-fit values of v and fl x o for both sets, focus- 
ing on the existence of some tighter constraints coming 
from the SNe la data. For the Hubble data we obtained 
a = 1.19 and £l x0 = 0.61 with x\oi = 0-786 whereas 
the SNe la data lead to v = 1.5 and tt x o = 0.70 with 
Xdof = 0.812 < 1; in broad terms the tighter constraints 
seems to be found with the Hubble data. Of course, these 
results can vary according to the parameter regions taken 
into account in the minimization process. 
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FIG. 1: The confidence contours associated with the la, 2a 
and 3a error bars are shown in the plane (v, Q x o)- The point 
[y, il x o) = (1.19, 0.61) shows the best-fit for the H ^ s (z) data. 



(a, p) 


(Qox ±a,v ±a) 


Xdof 


(1.01,0.05) 


(0.58 ±0.20, 1.19 ±1.13) 


0.786 


(1.01,0.1) 


(0.58 ±0.23, 1.19 ±1.13) 


0.786 


(1.01, 0.15) 


(0.61 ±0.02, 1.19 ±1.13) 


0.786 


(1.2,-0.1) 


(0.55 ±0.16, 1.19 ±1.13) 


0.786 


(1.2,-0.05) 


(0.57 ±0.17, 1.19 ±1.13) 


0.786 


(4/3,-0.1) 


(0.59 ±0.15, 1.19 ±1.13) 


0.786 


(4/3,0.1) 


(0.69 ±0.17, 1.19 ±1.13) 


0.786 



Zi 


H bs(zi) 


la 


redshift 


km s^ 1 Mpc^ 1 


uncertainty 


0.00 


74.2 


±3.6 


0.10 


69 


±12 


0.17 


83 


±8 


0.27 


77 


±14 


0.40 


95 


±17 


0.48 


97 


±60 


0.88 


90 


±40 


0.90 


117 


±23 


1.30 


168 


±17 


1.43 


177 


±18 


1.53 


140 


±14 


1.75 


202 


±40 



TABLE I: The observational H(z) data [g. 

Now, using the best-fit model parameters 6 C = 
(u,Q x0 ) = (1.19 ± 1.13,0.61 ± 0.02) we would like to 
compare the model having a MHRDE with the standard 
A-CDM scheme composed of baryonic matter and a con- 
stant dark energy flo x = 0.73 ± 0.04. In Fig.Q we show 
the H bs(z) data with their error bars, and the theoret- 



TABLE II: We show the observational bounds for the pair 
(Q. x o,u) varying (a,/3) with a given prior of Hq = 72.2 ± 
3.6 km s^Mpc -1 



ical curves corresponding to the A-CDM and MHRDE 
models, respectively. At small redshift < z < 0.2 both 
models seem to be in excellent agreement with the ob- 
servational H(z) data, while at large redshift z > 0.2 
the model shows a slight difference with the A-CDM. 
However, the observational data show that in the inter- 
val z £ [0, 1] both models have a good adjustment. As 
— 1 < uj(z) 1 uj x (z) < 0, the equations of state of the effec- 
tive fluid and dark energy do not cross the phantom line, 
at least for the best-fit model parameters used previously 
[see Fig.(@|]. The expression of uj x at present 



WzO 



v(a - /3)B 



(a-l)(l + v) + vB , 



(21) 



becomes u x o = —0.88 when evaluating at the best-fit 
values 6 C = (v,fl x0 ,a,l3) = (1.19,0.61,1.01,0.15). It is 
close to the value reported by WMAP-7, uj x o = —0.93, 
when the joint analysis of the WMAP+BAO±# ±SN 
data [3l[ for constraining the present-day value of the 
equation of state for dark energy is made. 
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FIG. 2: We show the confidence contour in the Ho — fi^o 
plane for the MHRDE model considering as priors v = 1.19, 
a — 1.01 and /3 — 0.15. The best fit value of (Ho, il x o) is used 
to compare with the ACDM model. 
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FIG. 3: We plot H(z) for the ACDM with H = 73.60 km/s 
Mpc and MHRDE model with prior H = 72.22 km/s Mpc. 
We also show the observational data H(zk) with their error 
bars. 



Fig. ([5]) shows the evolution of the decelerating param- 
eter q — —a/aH 2 with the redshift z for the MHRDE and 
A-CDM models. It takes the form 



'in 



l + i/(l-3B) 
2(1 + 1/) : 



(22) 



at z = for the former model. Using the best-fit values 
6 C = {v,n x Q,a,P) = (1.19,0.61,1.01,0.15) in Eq. ([2"2]). 
one gets qo = —0.27 while for the A-CDM model one 
obtains qo = —0.59. The critical redshift where the ac- 
celeration starts, 



1 + 



(2v- 1)5 



(l + i/)(l-B) 



1/3!/ 



(23) 



FIG. 4: We show the equations of state for the effective fluid 
u(z) and the dark energy ui x (z). 



turns to be z acc = 1.06 for the best-fit values 6 C , then 
our model enters the accelerated regime earlier than the 
ACDM one with z acc = 0.75. 

In Fig.© we plot the density parameters £l c , il x , its 
ratio r(z) and find the present-day values of fl x0 = 0.61, 
f2 C Q = 0.39 and r = 0.62. It shows that the model with a 
MHRDE seems to be appropriated for resolving the coin- 
cidence problem. Regarding the modified Ricci coupling 
function, one can show that Qm < and the coupling 
decreases its strength with the redshift and goes to zero 
in the far future, z — > — 1 [see Fig.©]. 

So far, we have discussed some cosmological con- 
straints coming from the observational H(z) data at dif- 
ferent redshifts. Let us take into account an alternative 
method, which is beyond parameter constraints, that can 
be implemented to distinguish the A-CDM from other 
dark energy models. The new diagnostic of dark energy 
called Om(z) is an interesting one because it does not 
involve the cosmic equation of state [33||. The Om(z) 
function is constructed from the Hubble parameter H 
and defined as [33 



Om(z) = 



H 2 



HI 



H*[(z + ir-lY 



(24) 



Since the Om(z) only depends upon the scale factor and 
its first derivative, one considers the present diagnostic 
method as a geometrical one. Futhcrmore, for the A- 
CDM model one finds that Om(z) = il m o so only if 
the dark energy is given by a cosmological constant the 
Om(z) remains constant. For our model we find that 
Om(z) can be written as 



Om(z) 



{B(l + zf + (1 - B){\ + z) 3 ^+^Y +1 - 1 



(Z + 1)3-1 



(25) 

Now, we are going to analyze some important regimes. 
In the far future, z — > —1, Eq. (j2"5"j) reaches the value 
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FIG. 5: The deceleration parameter q(z). 



FIG. 7: The Om(z) diagnostic. 
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FIG. 6: The density parameters (Q x , fl c ), the ratio r 
Cl c /£lx, and Hq Qm are shown versus the redshift z. 



1. For the present time, z — >• 0, the Eq. ([25]) gives 
[B + (1 - i?)(^ + l)]/(^ + 1) which is positive for */ > 0. 
In addition, using the constrained values of (y, f2 x o) = 
(1.19,0.61) coming from the H b s data, we plot the evo- 
lution of Om(z) for both the MHRDE and the A-CDM 
models [see Fig.©]. Interestingly enough, in the inter- 
val z G [0,1) the Om mhrde > Slmo the inequality is 
reversed at redshift, z > 1, as it happens with a phan- 
tom dark energy model. Thus the same behavior occurs 
with a quintessence model |33j . 

As a closing comment, we would like to address a dis- 
cussion concerning the values of a and /3 taken into ac- 
count through this section. Here we have focused on the 
transition of the Universe between a stage dominated by 
dark matter followed by an era dominated by the holo- 
graphic dark energy that makes the Universe exhibt an 
accelerated expansion (present-day scenario) and, in both 



stages a nonlinear interaction in the dark sector has been 
taken into account. In order to estimate the parameter 
v and Q x0 we have used the values of a and /3 which are 
consistent with the ^-statistical analysis because they 
fulfill the condition Xdof < 1- Now, we are going to ex- 
plore a modification on the aforesaid model by adding a 
third component, say p m , which does not interact with 
p c and p x . The total energy density reads as pt = p m + P 
with p = p c + Px and the MCE is split as 



p' + ap - (a - /3)p x 



0. 



ap„ 







(26) 



(27) 



so (|26|) is equal J2J while (|27)) shows that the third com- 
ponent is decoupled from the interacting dark sector, 
indicating that p m feels the presence of the dark sec- 
tor through the gravitational expansion of the Universe. 
From ([27]) one finds that p m = p m oa~ 3a whereas the 
behavior of dark matter and dark energy with the scale 
factor can be obtained from (fT4"]) and (TTol) . respectively. 
Using (Tl4"l) . (Tro|) and (|2"T|) one gets the Hubble parame- 
ter in term of the redshift x = 1 + z and the relevant 
cosmological parameters 



H(z) 

Hn 



(i-a m0 )(Bx 3 +(i-B)x 3 ^+V)^+n m0 x 3a 

(28) 



1/2 



m = v -^r L , i- n i- n „ n + (i - «)] ( 29 ) 

where the flatness condition 1 = Vl x q + £I c q + f2 m o has 
been used. In what follows we would to examine two 
traits of the model. First, we fix a = 4/3 to get a 
radiation contribution in the total density because it 
will address the problem of the dark energy at early 
times; thus, as is well known the fraction of dark energy 



7 



in the radiation era should fulfill the stringent bound 
fl x (z ~ 1100) < 0.1 in order for the model be consis- 
tent with the big bang nucleosynthesis (BBN) data. For 
the priors (H = 72.2, j; = 1.19, fl x0 = 0.7, a = 4/3) 
the Hubble data give as the best-fit values /? = 0.1 and 
O,„o = 1.7 x 10 6 along with x\ Q f = 0-78 < 1 with a frac- 
tion of dark energy £l x (z ~ 1100) = 0.2 nearly close to 
the BBN's bound. Second, employing the Hubble data 
for (f28|) we estimate the best-fit value of fi m o and a. Tak- 
ing as priors (H = 72.2, v = 1.19, Q x0 = 0.61, j3 = 0.15) 
the x 2_anar ysis yields as the best-fit values a = 1.01 
and fl m0 = 9.9 x 10~ 5 together with a Xdof = °-79 < 1 - 
Moreover, the latter case leads to an early dark en- 
ergy fl x (z ~ 1100) = 0.01 < 0.1 which is consistent 
with the bounds reported in [34J or with the future con- 
straints achievable by Planck and CMBPol experiments 
J35|. Therefore, taking the third component as the radi- 
ation term or a nearly radiation contribution, has helped 
to validate the first model, indicating that the value of 
the cosmological parameters selected are consistent with 
BBN constraints. 



IV. CONCLUSION 

We have investigated the interacting dark sector with a 
MHRDE, where the IR cutoff is provided by the modified 
Ricci scalar in a form of a linear combination of H and 
H 2 , so that p x = 2(H+3aH 2 /2)/(a-(3). This appealing 
cutoff yields a MCE in the dark sector with constant 
barotropic indexes, identified with a and j3. We have 
introduced an interaction between the dark matter and 
dark energy densities, homogeneous of degree 1 in the 
variables p and p', and solved the source equation for the 
total energy density of the mix. The equation of state 
of this effective fluid is that of the relaxed Chaplygin 
gas. The latter interpolates between a matter dominated 
phase at early times and an accelerated expanding phase 
dominated by the MHRDE at late times. 

We have used the observational Hubble data to con- 
strain the cosmological parameters of the model and to 
compare with the A-CDM model. Taking as a refer- 
ence point (af,/3f) — (1.01,0.15) we get the best fit at 
9 C = (v, Q x0 ) = (1.19, 0.61) with Xmi n = 7.86 leading to a 
good fit with Xdof = 0-786 < 1 per degree of freedom [see 
Fig- ([I])]. We have established that a model with a holo- 
graphic dark energy p x ex R leads to 0.59 < il x o < 0.69 
which is close to the bounds ft X Q = 0.73 provided by 
WMAP-7 [3l|. For (3 = 0.1, we have shown that the 



Ricci cutoff (a = 4/3) is consistent with other values of 
a because it fulfills the goodness condition {Xdof < -0- nl 
addition, wc have obtained the allowed range of (v, £l x o) 
when one varies a and /3 [see Table . pit ]. Properly esti- 
mating the Hq and £I x q with the Hubble data we have 
confronted the A-CDM with the MHRDE model; thus, 
both models give some bounds of the pair (Hq, £I x q) con- 
sistent the those reported in [31|. Besides, we have taken 
into account the SNe la with the Union2 data for calcu- 
lating the best-fit values of v and fl x o. It led to v = 1.5 
and n x0 = 0.70 with x 2 dof = 0.812 < 1 while the Hubble 
data gave v = 1.19 and fl x o = 0.61 with Xdof = 0.786. 

Wc have found that the equations of state of the dark 
energy equation and the unified fluid, at the best-fit val- 
ues 9 C , do not cross the phantom divide line [see Fig.Q] 
while the present value of the equation of state for the 
dark energy is lu x q = —0.88. 

From the deceleration parameter [see Fig. ([5])] and the 
best fit values 8 C , we have obtained that the acceleration 
starts at z acc = 1.06 hence, the model with a MHRDE 
enters the accelerated regime earlier than the A-CDM 
with z^c = 0.75. We have shown that the density pa- 
rameters il c , tt x , and its ratio r(z) in Fig.© seem to 
alleviate the coincidence problem. It is related to the 
decreasing behavior of the interaction with the redshift 
and its vanishing limit in the far future [see Fig.©]. Wc 
have applied the geometrical Om diagnostic method to 
our model which seems include several aspects different 
from the A-CDM model. 

Finally, we have analyzed a new model with a third 
fluid decoupled from the dark interacting sector. The 
latter scheme allowed us to examine the issue of the early 
dark energy at z ~ 1100 as well as to check the match of 
the models with the BBN constraints [3J], |35| . 
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